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Abstract 


Turbulent convection is an inherently non-local phenomenon and a primary condition 
for a successful treatment of the CBL (convective boundary layer) is a reliable model of 
non-locality. In the dynamic equations governing the convective flux, the turbulent kinetic 
energy, etc., non— locality is represented by the third— order moments, TOMs. Since the 
simplest form, the so-called down gradient approximation (DGA), severely underestimates 
the TOMs (up to an order of magnitude), a more physical model is needed. In 1994, an 
analytical model was presented which was derived directly 'rom the dynamical equations 
for the TOMs. It considerably improved the DGA but was a bit cumbersome to use and, 
more importantly, it was based on the quasi— normal iQN) approximation for the 
fourth-order moments. 

Here, we present a new analytic expression for the TOMt which is structurally simpler 
than the 1994 expression and which avoids the QN approximation. The resulting fit to the 
LES data is superior to that of the 1994 model. 
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I. Introduction. 

The search for a reliable expression for the third-order moments to be used in the 
dynamic equations for the second-order moments such as the turbulent kinetic energy, the 
convective fluxes, etc., has a long history. For many years, people used the so-called down 
gradient approximation but LES (large eddy simulations) have shown that said model 
severely underestimates the TOMs (Moeng and Wyngaard, 1989). 

Prompted by these results, Canuto et al.(1994) undertook the task of solving directly 
the dynamic equations for the TOMs thus avoiding the need for phenomenological 
expressions. The key merit of the 1994 model was to exhibit the fact that all the TOMs are 
a linear combination of the gradients of all the second— order moment j and not only of 
selected ones, as assumed in the down-gradient approximation. Fron the performance 
viewpoint, the new TOMs reproduced the LES data quite satisfactorily but the predicted 
w 2 !? and w# 2 were not as good as that of the other TOMs. The weakest point in the 1994 
model was the use of the QN (quasi-normal) approximation for the fourt l-order moments. 
Both these limitations motivated us to search for new expressions for the TOMs which are 
simpler and with a better physical content. 

II. The new physical ingredient of the third— order moments 

There are six TOMs to be considered: 

w 3 , cpw, w 2 "#, w# 2 , ’5 3 , q 2 ^ (1) 

Here, u,v,w and 9 are the fluctuating velocity and temperature fields and q 2 =u 2 +v 2 +w 2 . 
Since the TOMs in (1) have different dimensions, we multiply the last four by appropriate 
variables so that all the TOMs have dimensions of a velocity cubed. Thus, we introduce the 
new variables x's which have the same dimensions: 


X 1 = ga^w 2 ^, * 2 = (g or v ) 2 w^ 

(2a) 

* 3 = (ga^) 3 # 3 , x^ 

(2b) 

X = wq 2 , Z E w 3 

(2c) 
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Here, a is the volume expansion coefficient, g is the local gravity and a time scale that 
will be discussed below. The original dynamic equations for the TOM given in (1) can be 
found in Canuto (1992), Eqs. (37a), (38a), (39a) and (40a). These equations entail fourth 
order moments which can be written in general as: 

abed = (ab cd + ac M 4- acf 5c)f (3a) 

When the function f is taken to be unity, Eq.(3a) corresponds to the quasi-normal (QN) 
approximation. 'Jnder the QN approximation, the time scale introduced before can only 
be identified with the dynamic time scale of turbulence r= 2K / e (K and c are the turbulent 
kinetic energy and its rate of dissipation). 

With r =r, the explicit form of the TOMs was presented in Canuto et al.( ? 994). 

V 

Even though the comparison with LES data was overall satisfactory, the predicted w 2 ? and 
w^ (figs. 10—11 of Canuto et ah, 1994) were not as good as that of the other TOMs. In 
addition, the expressions of the TOMs were rather cumbersome to use due to the fully 
explicit form of the determinant obtained when solving the algebraic set of equations for 
the TOMs. 

For these reasons, we felt motivated to search for new expressions for the TOMs that 
are both simpler to handle and possess a better physical content. Zilitinkevich et al. (1999) 
tried to do so but their solution is not practical since it does not provide an expression for 
w 3 which must be derived form outside, say from an LES. 

To improve the 1994 model, we were guided by the fact that with f=l, the TOM can 
become arbitrarily large while in reality they have finite values, e.g., the small value of the 
skewness. Thus, reducing the TOM given by the f=l case is a way to avoid unphysical 
results. How to directly relate this "damping effect" to a f^l is not a matter that can be 
carried out analytically, rather, we have used f in (3a) to formally indicate the physical 
motivation of our new approach. In the most successful heuristic model used to cut down 
the growth of the TOM, the EDQNM model (Lesieur, 1992), the damping is represented by 
an additional time scale which one must choose on physical grounds. Here, we suggest to 
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t = rfl+A NV] 4 , N 2 =-gadT/dz, r=2K/e 
v o 

We shall further take: 

A =0.04 if N 2 >0, A =0 if N 2 <0 
o o 

The previous model (Canuto et ah, 1994) corresponds to f=l and/or to A o =0. 

III. The new expression for t >e third-order moments 

The new analytic expressions for the TOMs are quite simple: 

x = X z-X 

1 o 1 

X = Y z- Y 

2 0 1 

x = Z x — Z 

3 0 2 1 

x = W x + c 4 x + W 

4 0 5 2 1 

x = 9 z-fi 
5 0 1 

z = (n - i.2x - k )(c - i.2x + n y 1 

V 1 1 2 5 /v o o' 

where the functions X, Y, Z, W and f l are defined as: 

x 0 = 7 2 ^ 2 (l-7 3 ^ 2 )[l-(7-j-7 3 )N 2 ]' 1 

X = [7 f + 7 f + 7 (1-7 ft 2 )f ][1 — (7+7 )N 2 ]' 1 
1 1 'o 0 11 2 K 3 ’ 2 Jl 1 3' 

Y = 27 ft 2 (l-7 ft 2 ) 4 X 

0 2 3 0 

Y = 27 (1-7 ft 2 ) 4 (ft 2 X + 7 7 _1 f + f ) 

1 V 3 ; 1 0 1 0 1 

Z = |(c-2) 4 ft 2 
Z = |(c— 2)' 1 f 

1 2 V 0 

Q = u X + o;Y 
0 00 10 

fl = til X + Y + LJ 
1 0 1 11 2 

W = i ft 2 , w = - c 4 f 
0 2 C 1 3 

The auxiliary functions u/s are: 

u = 7 (1— 7 ft 2 )' 1 , oj = (2c) 4 ui , u) — u> f + j u) f 
Finally, the 7's are constants which depend on the only adjustable parameter 



( 7 ) 


7 = 0.52c' 2 (c-2)' 1 , 7 = 0.87c‘ 2 7 =0.5c' 1 

'0 v ' 1 2 

7 = otoc - 1 ^)- 1 
3 

7 = 2.4(3c+5) 4 , 7 = O.Gc-^c+b)- 1 

4 5 

Based on previous work, the suggested value is c=7 but small variations are allowed. The 


second order moments enter through the functions f which are defined as follows: 

0v>5 

r / \ 3 4 T 

f „ = V8“) r v J ® 

f , = <« a ) 2r J( J § + > 


f 2 = 8< J W + 2 S aT l**% 

r 2 / — 2<5J , t <9K \ 

f = gar 2 (w 2 ^ + J w ) 


dJ 


f = r^( 

. \r ■ 


■ dvP , <9K 


) 


f = r W 2 


0W 2 


5 V 

All the functions fs have dimensions of velocity cubed. Finally: 


J=W3, N 2 = r 2 N 2 


( 8 ) 

(9) 


IV. Test of the new TOM vs LES data. 

In Fig.l we compare the new TOMs given by Eqs.(2) and (4) vs. LES data. Rather 
than solving the CBL dynamic equations as it was done in the 1994 paper, here we 
employed LES data to compute the second— order moments Eq.(8), r=2K/e and N 2 (z). 

The better agreement with LES data with respect to the 1994 model, especially in 
Figs lc-d, is also partly due to the simpler analytical form of the TOM which has allowed 
us to test slight variations around the c=7 value. Due to its rather rigid nature, the 1994 
model did not allow the same freedom. However, the key reason for the better performance 
of the new model is of physical origin: we have abandoned the quasi— normal approximation 
for the fourth-order moments that we employed in the 1994 model. This means that we 
have searched for a way to cut down an otherwise unphysical growth of the TOM by 
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adopting an EDQNM-like procedure. Regrettably, at present we dont have an a priori 
derivation for (3b, c) which must thus be considered an heuristic suggestion. 

V. Conclusions. 

The results presented in Fig.l satisfy the two requirements set out at the beginning, 
the expressions for the TOMs are simpler than those of the 1994 model and their physical 
content is better. As a result, the large values of w 2 ? and w# 2 that characterized the 1.994 
model are no longer present and an overall better fit is obtained. 
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Figure Caption 

Fig.l. The third— order moments Eq.(2) vs. z/h. The LES results of Zilitinkevich et al. 
(1999) are plotted as dotted lines (the LES data did not contain the q 2 '# value); the 
down— gradient approximation (DGA) model is plotted as dashed lines while the present 
model results are plotted as solid lines. As well known, the DGA severely underestimates 
the third-order moments. All TOMs are normalized with Deardorff convective scales 
w*(=2ms 4 ) and 0*(=O.12K). The value of the PBL depth h is 10 0m. 
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